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Abstract 



Microscopic black holes explode with their temperature varying inversely 
as their mass. Such explosions would lead to the highest temperatures in 
the present universe, all the way to the Planck energy. Whether or not a 
quasi-stationary shell of interacting matter undergoing radial hydrodynamic 
expansion surrounds such black holes is controversial. In this paper relativistic 
viscous fluid equations are applied to the problem assuming sufficient particle 
interaction. It is shown that a self-consistent picture emerges of a fluid just 
marginally kept in local thermal equilibrium; viscosity is a crucial element of 
the dynamics. 
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Hawking radiation from black holes [|l| is of fundamental interest because it relies on the 
application of relativistic quantum field theory in the presence of the strong field limit of 
gravity, a so far unique situation. It is also of great interest because of the temperatures 
involved. A black hole with mass M radiates thermally with a temperature 



where mp = G"^/^ = 1.22 x 10^^ GeV is the Planck mass. (Units are h = c = k-Q = 1.) In 
order for the black hole to evaporate it must have a temperature greater than that of the 
present-day black-body radiation of the universe of 2.7 K = 2.3x10^"^ eV. This implies that 
M must be less than 1% of the mass of the Earth, hence the black hole most likely would 
have been formed primordially and not from stellar collapse. The black hole temperature 
eventually goes to infinity as its mass goes to zero, although once Th becomes comparable to 
the Planck mass the semi-classical calculation breaks down and the regime of full quantum 
gravity is entered. Only in two other situations are such enormous temperatures achievable: 
in the early universe (T similarly asymptotically high) and in central collisions of heavy nuclei 
like gold or lead (T = 500 MeV is expected at the RHIC (Relativistic Heavy Ion Collider) 
just completed at Brookhaven National Laboratory and T = 1 GeV is expected at the LHC 
(Large Hadron Collider) at CERN to be completed in 2005). The spontaneously broken 
chiral symmetry of QCD gets restored in a phase transition/rapid crossover at a temperature 
around 160 MeV, while the spontaneously broken gauge symmetry in the electroweak sector 
of the standard model gets restored in a phase transition/rapid crossover at a temperature 
around 100 GeV. The fact that temperatures of the latter order of magnitude will never be 
achieved in a terrestrial experiment motivates me here to study the fate of primordial black 
holes during the final minutes of their lives when their temperatures have risen to 100 GeV 
and above. The fact that primordial black holes have not yet been observed does not 
deter me in the least. 

There has been some controversy over whether the particles scatter from each other after 
being emitted, perhaps even enough to allow a fluid description of the wind coming from 
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the black hole. First the situation will be reviewed, and then it will be shown that a self- 
consistent description is afforded by the application of relativistic viscous fluid equations, at 
least towards the end of the black hole's existence and under certain assumptions. 

When Th -C me (electron mass) only photons, gravitons, and neutrinos will be created 
with any significant probability. These particles will not interact with each other but will be 
emitted into the surrounding space with the speed of light. Even when Th ~ rUe the Thomson 
cross section is too small to allow the photons to scatter very frequently in the rarified 
electron-positron plasma around the black hole. This may change when Th ~ 100 MeV 
when muons and charged pions are created in abundance. At somewhat higher temperatures 
hadrons are copiously produced and local thermal equilibrium may be achieved, although 
exactly how is an unsettled issue. Are hadrons emitted directly by the black hole? If so, 
they will be quite abundant at temperatures of order 150 MeV because their mass spectrum 
rises exponentially (Hagedorn growth as seen in the Particle Data Tables ||^). Because they 
are so massive they move nonrelativistically and may form a very dense equilibrated gas 
around the black hole. But hadrons are composites of quarks and gluons, so perhaps quarks 
and gluon jets are emitted instead? These jets must decay into the observable hadrons on 
a typical length scale of 1 fm and a typical time scale of 1 fm/c. Once the hadrons appear 
they may form an equilibrated gas around the black hole just as if they had been produced 
directly albeit with some time delay. One can find arguments both for and against 
1^ thermal equilibrium among the strongly interacting hadrons outside the Schwarzschild 
radius. Recently a numerical study of photosphere formation around a primordial black 
based on a relaxation time approximation to the Boltzmann equation has been performed 
for QED and QCD interactions 0. It was found that significant particle scattering would 
lead to a photosphere though not perfect fluid flow. Certainly this is a very difficult and open 
problem in quantum statistical mechanics, just as it is in high energy heavy ion collisions 

i- 

Let us assume that a primordial black hole is surrounded by a shell of expanding in- 
teracting matter in approximately local thermal equilibrium when Th is large enough. A 
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detailed description of how this situation comes to be is a difficult problem as discussed above 
and is not addressed in this paper. The relativistic imperfect fluid equations describing a 
steady-state, spherically symmetric flow with no net baryon number or electric charge and 
neglecting gravity (see below) are T^'^.^^ = black hole source. The nonvanishing components 
of the energy-momentum tensor in radial coordinates are M 



diss 

T- = t;V(i^ + e)+i^ + ATdiss (2) 

representing energy density, radial energy flux, and radial momentum flux, respectively, in 
the rest frame of the black hole. Here v is the radial velocity with 7 the corresponding 
Lorentz factor, u = vf, e and P are the local energy density and pressure, and 

ATd- = -- 72 f — - - f — + —\ (3) 

''''' 3 \dr r J \dr r J ^ 

where rj is the shear viscosity and ( is the bulk viscosity. A thermodynamic identity gives 
Ts = P + e for zero chemical potentials, where T is temperature and s is entropy density. 
There are two independent differential equations of motion to solve for the functions T(r) 
and v{r). 

An integral form of these equations is probably more useful since it can readily incor- 
porate the input luminosity Lj from the black hole. The flrst represents the equality of the 
energy flux passing through a sphere of radius r with the luminosity of the black hole. 

ATTr^T^"- = Li (4) 

The second follows from integrating a linear combination of the differential equations. It 
represents the combined effects of the entropy from the black hole together with the increase 
of entropy due to viscosity. 



The term Li/Th arises from equating the entropy per unit time lost by the black hole —dS\i/dt 
with that flowing into the matter. Using the area formula for the entropy of a black hole, 
'S'h = mpTTTg = iirM"^ /rup, and identifying —dM/dt with the luminosity, the entropy input 
from the black hole is obtained. 

The above pair of equations are to be applied beginning at some radius rj greater than 
the Schwarzschild radius r^, that is, outside the quantum particle production region of the 
black hole. The radius at which the imperfect fluid equations are flrst applied should be 
chosen to be greater than the Schwarzschild radius, otherwise the computation of particle 
creation by the black hole would be invalid. It should not be too much greater, otherwise 
particle collisions would create more entropy than is accounted for by eq. (5). The energy 
and entropy flux into the fluid come from quantum particle creation by the black hole at 
temperature T/j, and particle production is dominated by particles with mass less than the 
temperature. Massless particles emitted from a surface at rest have an average outward 
velocity of l/VS. Therefore it is not unreasonable to assume the initial flow velocity is 
Vi = v{ri) = l/\/3. Furthermore, considerations of smoothness and continuity suggest that 
dv/dr — at r^. In fact one should not expect the precise choice of initial conditions to 
matter at large radii; see the discussion of the scaling solutions later in this paper. Once 
the functions s(T), "^(T), and C(^) ^"^^ specified, and Tj = T{ri) can be determined by 
the integral form of the equations of motion, eqs. (4) and (5). Gravitational effects are of 
order rs/r, hence negligible for r > (5 — 10)r5. 

A black hole has a Schwarzschild radius rs = 2M/m| = l/AnTh. Note that nTh • 2rs = 
1/2. Roughly, the average thermal momentum of a massless particle times the diameter of 
the black hole is 1/2. This is just a manifestation of the uncertainty principle applied to the 
creation of an excitation in a confined region of space. The luminosity is 

L = -^ = a{M)^^=64nMTM (6) 

where a{M) is a function reflecting the species of particles available for creation in the grav- 
itational field of the black hole. It is generally sufficient to consider only those particles with 
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mass less than T^, more massive particles are exponentially suppressed by the Boltzmann 
factor. Then 

a = 2.011 X 10"® [4200Aro + 2035Ari/2 + SSSA^i + gSA^a] . (7) 

Here A^^ is the net number of polarization degrees of freedom for all particles with spin 
s and with mass less than T^. The coefficients for spin 1/2, 1 and 2 were computed in 
ref. 1^ and for spin in ref. [|lOl- In the standard model Nq = 4 (Higgs), N1/2 = 90 
(three generations of quarks and leptons), Ni = 24 (SU(3)xSU(2)xU(l) gauge theory), and 
N2 = 2 (gravitons). This assumes Th is greater than the temperature for the electroweak 
gauge symmetry restoration [jlT|. Numerically a{Th > 100 GeV) = 4.43 x 10^^. Starting 



with a black hole of temperature T^, the time it takes to evaporate/explode is 

rup , . 

^=3a(T,)(87rT,)3- 

This is also the characteristic time scale for the rate of change of the luminosity of a black hole 
with temperature T^. Generally At ^ l/7)i, thus justifying the quasi- stationary assumption 
|T^. For example, a black hole with temperature 1 TeV has a Schwarzschild radius of 



1.57 X 10^^ fm, a mass of 10^° g, a luminosity of 7 x 10^'' erg/s, and has 464 seconds to 
live. The input luminosity to the expanding fluid Li will be less than the total luminosity L, 
because gravitons will escape without scattering, and neutrinos could scatter only in those 
regions where the temperature is greater than about 100 GeV. The reason is neutrino cross 
sections above 100 GeV will be similar to the cross sections of other fermions because of the 
restoration of the spontaneously broken electroweak gauge symmetry. 

Determination of the equation of state as well as the two viscosities for temperatures 
ranging from MeV to TeV and more is a formidable task. Here we shall consider two 
interesting limits and then a semi-realistic situation. A realistic, quantitative description of 
the relativistic black hole wind, including the asymptotic observed particle spectra, is left 
for a future pubhcation. 

First, consider the adiabatic limit (like milk) with an equation of state e = aT^, s = 
{A/3)aT^, and rj = ( = 0. This is equivalent to assuming that the mean free paths of the 



particles are all small compared to the length scale over which the temperature and other 
thermodynamic quantities change significantly. A scaling solution, valid when 7 ^ 1, is 
T(r) = To(ro/r) and 7(r) = 7o(r/ro), where 70T0 = Th- The tq is any reference radius 
satisfying the stated criterion. 

Second, consider the highly viscous isoergic limit (like honey) in the sense that the flow 
velocity approaches a limiting value vo at large r. This requires a power-like equation of state 
e (X and viscosities rj oc ( (x T^l'^ . It results in the scaling solution T(r) = To(ro/r)^/^. 
This is not very realistic: a massless gas with dimensionless coupling constants and 5 = 4 
would require viscosities of order whereas one would expect on dimensional grounds. 

Now consider a semi-realistic situation with e = aT^, s = (4/3)aT^, 77 = bsT^, and ( = 
hsT^- This is typical of relativistic gases with dimensionless coupling constants, although 
quantum effects will give logarithmic corrections |jl^,|13|]. A scaling solution, valid at large 
radii when 7 ^ 1, is T{r) = TQirQ/rY^^ and 7(r) = '^^{r /roY/"^ . The constants are related by 
SGaTgro = (3265 + 4416^)70. This r-dependence of T and 7 is exactly what was conjectured 
in ref . 0| . 

Is the semi-realistic situation described above really possible? Can approximate local 
thermal equilibrium, if achieved, be maintained? The requirement is that the inverse of the 
local volume expansion rate 9 = u^.^ be comparable to or greater than the relaxation time 
for thermal equilibrium P]. Expressed in terms of a local volume element V and proper 
time r it is ^ = {l/V)dV/dT, whereas in the rest frame of the black hole the same quantity 
can be expressed as {l/r'^)d{r'^u)/dr. Explicitly 



Of prime importance in achieving and maintaining local thermal equilibrium in a relativistic 
plasma are multi-body processes such as 2 — 3 and 3 — > 2, etc. This has been well-known 
when calculating quark-gluon plasma formation and evolution in high energy heavy ion 



This is a formidable task in the standard model with its 16 species of particles. Instead 
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(9) 



collisions []T3[ and has been emphasized in ref. in the context of black hole evaporation. 



we make three estimates for the requirement that local thermal equilibrium be maintained. 
The first and simplest estimate is to require that the thermal DeBroglie wavelength of a 
massless particle, 1/3T, be less than 1/9. The second estimate is to require that the Debye 
screening length for each of the gauge groups in the standard model be less than 1/9. The 
Debye screening length is the inverse of the Debye screening mass where n = 1, 2, 3 
for the gauge groups U(l), SU(2), SU(3). Generically m° oc gnT where gn is the gauge 
coupling constant and the coefficient of proportionality is essentially the square root of the 



number of charge carriers [T^. For example, for color SU(3) = g-^Jl + Nf/6T where 



A^f is the number of light quark flavors at the temperature T. The numerical values of the 
gauge couplings are: gi = 0.344, g2 = 0.637, and g^ = 1.18 (evaluated at the scale mz) 0. 
So within a factor of about 2 we have mP ^ T. The third and most relevant estimate is 
the mean time between two-body collisions in the standard model for temperatures greater 
than the electroweak symmetry restoration temperature. This mean time was calculated 



in ref. [|T^ in the process of calculating the viscosity in the relaxation time approximation. 
Averaged over all particle species in the standard model one may infer from that paper an 
average time of 3.7/T. Taking into account multi-body reactions would decrease that by 
about a factor of two to four. All three of these estimates are consistent within a factor of 
2 or 3. The conclusion to be drawn is that local thermal equilibrium should be achieved 
when 9 < T, or 84a/{32bs + 4416^) < 1. This criterion places a minimum value on the 
strength of the viscosities. Once thermal equilibrium is achieved it is not lost because 9/T 
is independent of r. The picture that emerges is that of an imperfect fluid just marginally 
kept in local equilibrium by viscous forces. 

The hot shell of matter surrounding a primordial black hole provides a theoretical testing 
ground rivaled only by the big bang itself. In addition to the questions already raised, one 
may contemplate baryon number violation at high temperature and how physics beyond 
the standard model might be important in the last few minutes in the life of a primordial 
black hole. Finally, such black holes may contribute to the highest energy cosmic rays whose 
origin is a long-standing puzzle. Experimental discovery of exploding black holes will be one 
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of the great challenges in the new millennium. 
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